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Introduction 



Consider a system of n competing species whose states are characterized by vectors in 
the closed positive cone K := [0, oo) n C R n . When time is discrete the development of 
the system is given by a continuous map T: K — > K. For continuous time is continuous 
the development is governed by a periodic system of differential equations x = F(t, x) = 
F(t + l,x). In this case T denotes the Poincare map. 

For discrete time the trajectory of a state x is the sequence {T k x}, also denoted by 
. {x(k)}, where k varies over the set N of nonnegative integers. In the case of an autonomous 

differential equation (i.e., F is independent of t), the trajectory of x is the solution curve 
through x, denoted by T l x or x(t), where t 6 [0, oo). In both cases the limit set lo(x) is the 
set of limit points of sequences x(tfc) where t& — ► oo. 

In order to exclude spontaneous generation we assume Tj(x) = when X{ = 0. Thus 
there are functions Gi : K — > [0, oo), assumed continuous, such that 

Ti(x) = XiGi(x), (x € K, i = 1,... ,n) (1) 

For continuous time we assume the differential equation is a system of having the form 
Xi = XiGi(t,x). If Xi is interpreted as the size of species i then Gi(x) is its per capita 
growth rate. 

We take "competition" to mean that increasing any one species does not tend to increase 
the per capita growth rate of any other species, conventionally modeled by the assumption 

§t<°> 

A carrying simplex for the map T is a compact invariant hypersurface S C K such 
that every trajectory except the the origin is asymptotic with a trajectory in E, and S is 
unordered for the standard vector order in K. In the case of an autonomous differential 
equation we require that S be invariant under the maps T* for all t > 0. Some maps 
have no carrying simplices, others have infinitely many. Our main results gives conditions 
guaranteeing a unique carrying simplex. 
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Terminology 

A set Y C K is positively invariant under a map or an autonomous differential equation if 
it contains the trajectories of all its points, so that T l Y C Y for all t > (here t € N or 
[0, oo) as appropriate). We call Y invariant if T l Y = Y for all t > 0. 

If S is a differentiable map, its matrix of partial derivatives matrix at p is denoted by 
S'(p). 

The geometry of K plays an important role. For each subset I C {1, . . . , n} the I'th 
facet of K is 

Ki = {x G K: X j = j i 1} 

Thus K{j} is the i'th positive coordinate axis. A facet is proper if it lies in the boundary of 
K, meaning I / {0}. The intersection of facets is a facet: Kj n K j = K/yj. The boundary 
of K in R n , denoted by K, is the union of the proper facets. Each x € K\{0} belongs to the 
unique facet Kj^) where I(x) := {i: Xi> 0}. 

For each nxn matrix A and nonempty I C {1, . . . , n} we define the principal submatrix 

A i ■= IMjei 

The vector order in R n is the relation is defined by x >z y x — y € K. We write 

x >- y if also x ^ y. For each set I C {1, ... , n} we write x >~i y if x, y G Ki and x >- y, and 
x y if also Xj > yi for all igl. The reverse relations are denoted by ^, -< and so forth. 

The closed order interval defined by a, b € R n is 

[a, b] := {x E R ra : a ^ x r< 6} 

Carrying simplices 

A carrying simplex is a set S C K\{0} having the following properties: 
(CS1) S is compact and invariant. 

(CS2) for every x € K\{0} the trajectory of x is asymptotic with some y € S, i.e., 
lim^oo |T*x - T*y| = 0. 

(CS3) S is unordered: if x, y G S and x ^ y then x = y. 

It follows that each line in K through the origin meets S in a unique point. Therefore 
E is mapped homeomorphically onto the unit (n-l)-simplex 

A"- 1 :={xGK: J^x, = 1} 

i 

by the radial projection x i— > x/(^ Xj). 

Long-term dynamical properties of trajectories are accurately reflected by the dynamics 
in S by (CS1) and (CS2), and (CS3) means that S has simple topology and geometry. The 
existence of a carrying simplex has significant implications for limit sets lu(x): 
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• If x > then u(x) C X, a consequence of (CS2). In particular, S contains all nontrivial 
fixed points and periodic orbits. 

• If a, b € K are distinct limit points of respective states x, y >z (possibly the same 
state), then there exist i,j such that cij > 6j, Oj < this follows from (CS3). Thus 
either u{x) = uj(y), or else there exist i,j such that 

limsup Xi(t) — yi{t) > 0, liminf Xj(t) — yj(t) < 

In many cases S is the <?Zo6a/ attractor for the dynamics in K\{0}, meaning that as 
t goes to infinity, the distance from x(t) to £ goes to zero uniformly for x in any given 
compact subset of K\{0}. This implies (Wilson [35]) that there is a continuous function 
V: K\{0} -> [0, oo) such that if x / then 

• V"(x) = xGS, 

• V(x{t)) < V(x) x <£ E, 

• lim^oo K(x(t)) = 0, 

We can think of V as an "asymptotic conservation law". While there are many such 
functions for any carrying simplex, it is rarely possible to find a formula for any of them. 
Before stating results we give two simple examples for n = 1: 

Example 1 

If T is the time-one map for the flow defined by the logistic differential equation 

x = rx{a — x), r,a>0, {x > 0), 

the carrying simplex is just the classical carrying capacity a. Here one can define V{x) = 
\x — a\ for x > 0. 

Example 2 

Consider the map 

T: [0, oo) — >[0, oo), Tx = xe b ~ ax , b,a>0, x £ [0, oo) (2) 

Note that 

T'{x) = (1 - x)e b - ax , T'(b/a) = 1-6 

If there is a carrying simplex, it has to be the unique positive fixed point b/a, in which case 
Hindoo T k x = b/a for all x > 0. 

Ifb < 1 then b/a is the carrying simplex. In this case the maximum value of T is taken 
uniquely at 1/a > b/a. If < x < b/a then x < Tx < b/a, hence T k x — > 6/a. It follows 
that if some T J 'a; < b/a then again T fc x — > 6/a. If the entire orbit of x is > b/a then the 
sequence {T fc x} decreases to a fixed point > b/a, hence to be b/a. 

If b > 2 i/iere is no carrying simplex. For then |T'(6/a)| > 1, making b/a a locally 
repelling fixed point. The only way the trajectory of y ^ b/a can converge to b/a is for 
T^y = b/a for some j > 0. The set of such points y is nowhere dense because T is a 
nonconstant analytic function, hence there is no carrying simplex. For sufficiently large b 
the dynamics is chaotic. 

Example 5, below, is an n-dimensional generalization of Equation (2). 
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We say that T is strictly sublinear in a set X C K if the following holds: x G X and 
< A < 1 imply \x <E X and 

XT(x) -< T(Xx), (xeX\0) (3) 

Thus the restricted map T\X exhibits what economists call "decreasing returns to scale." 
A state x majorizes a state y\ix>-y, and x strictly majorizes y if Xi > implies Xi > yi. 
The map T: K — > K is strictly retrotone in a subset X C K if for all ijelwe have 

majorizes Ty ==> x strictly majorizes y 

Equivalently: 

i,t/6lnKi and Tx >- Ty ==> x 3>i y 

The origin is a repellor if T _1 (0) = and there exists 6 > and an open neighborhood 
II 7 C K of the origin such that liminf^oo \T k x\ > 5 uniformly in compact subsets of W\{0}. 

If in addition there is a global attractor T, as will be generally assumed, then T contains 
a global attractor Tq for T| K\{0}. In 

We will assume T is given Equation (1) has the following properties: 

(CO) T-^O) = and Gi(0) > 1. 

The first condition is means that no nontrivial population dies out in finite time. The 
second means that small populations increase. 

(CI) There is a global attractor T containing a neighborhood o/0. 

Together with (CO) this implies that there is a global attractor To C T for T| K\{0}. 
The connected component of the origin in K\ro is the repulsion basin B(0). 

(C2) T is strictly sublinear in a neighborhod ofT. 
This holds when < A < 1 G(x) -< G(Xx). 

(C3) T is strictly retrotone in a neighborhood of the global attractor 
A similar property was introduced by Smith [29]. 

Denote the set of boundary points of V in K by d^T. 

Theorem 3 When (CO) — (C3) hold, the unique carrying simplex is S = <9|<r = <9«-B(0), 
and S is the global attractor for T\ K \ {0}. 

The proof will appear elsewhere. 

The same hypotheses yield further information. It turns out that if T\T is locally 
injective (which Smith assumed), it is a homeomorphism of T; and in any case the following 
condition holds: 

(C4) The restriction of T to each positive coordinate axis K^j has a globally attracting 
fixed point . 
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We call c/(j) an axial fixed point. Denoting its i'th coordinate by qi > 0, we set 

q := (qi,...,q n ) = 

i 

Smith [29] shows that (C3) and (C4) imply (CI) with T C [0,q]. In many cases the easiest 
way to establish a global attractor is to compute the axial fixed points and apply Smith's 
result. 

The following condition implies (C3) for maps T having the form (1) when G is C 1 : 
(C5) If x £ K\{0}, the matrix [G'^x)]^^ has strictly negative entries 

For d G R n we denote the diagonal matrix D with diagonal entries Da := d{ by [<i] diag 
and also by [d;] diag . The n x n identity matrix is denoted by /. 
A computation shows that 

T'(x) = [G(x)f as + [xf ae G'{x). 

When x is such that all Gi{x) > 0, this can be written 

T'(x) = [G(x)] d ^(I-M(x)), 

x 



Mix) ■ LG.W 

and the entries in the n x n matrix M(x) are 



1 d'ag (4) 

G'(x), 



Gi(x) oxj 

d\ogG i ( . 
= -Xi—R {x) 

OX i 



(5) 



Note that (C5) implies My(x) > 0. 

The spectral radius p{M) of an n x n matrix M is the maximum of the norms of its 
eigenvalues. It is a standard result that if p{M) < 1 then I—M is invertible and (I— M) _1 = 

Theorem 4 Suppose G is C 1 . Assume (CO), (CI), (C2), (C5), let (C4) hold with [0,q] C 
X, and assume 

O^x^q => p{M(x)) < 1 (6) 

Then (C3) holds, whence the hypotheses and conclusions of Theorem 3 are valid. 

The proof will be given elsewhere. Under the same hypotheses the following conclusions 
also hold: 



T\T is a diffeomorphism 

i/iSTnKJ then the matrix [T'(x)i] _1 has strictly positive entries. 
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When (C5) holds, either of the following conditions implies (6): 

0<x±q Y^Mij(x)<l, (j = l,...,n) (7) 

i 

~> x <q ^2Mij(x) <1 (i = l,...,n) (8) 

Each of these conditions implies that the largest positive eigenvalue of M(x) is the spectral 
radius by (C5) and the theorem of Perron and Frobenius [2], and that this eigenvalue is 
bounded above by the maximal row sum and the maximal column sum by Gershgorin's 
theorem [3]. 



Competition models 

In the following illustrative examples we calculate bounds on parameters that make row 
sums of M{x) obey (7), validating the hypotheses and conclusion of Theorem 4 and 3. 

Example 5 

Consider a multidimensional version of Equation (2), based on an ecological model of May 
k Oster [22]: 

T: K — K, Ti(x) = x { exp (B { - ^ A ijXj ) , B h A {j > (9) 

This map is not locally injective. In a small neighborhood of the origin T is approximated 
by the discrete-time Lotka-Volterra map T defined by (Tx)i = (exp Bj)xi(l — ^ - A^Xj), 

but as T does not map K into itself, it is not useful as a global model. T has a global 
attractor V and a source at the origin, so a carrying simplex is plausible. But the special 
case n = 1, treated in Example 2, shows that further restrictions are needed. 

Condition (C5) holds with Gi(x) = exp — ^ - A^Xj^. Evidently these functions are 
strictly decreasing in x, which implies T is strictly sublinear. (C4) holds with = Bi/Aa, 
and it can be shown that V C [0, q\. In (4) the matrix entries are 

M ij (x)=x i A ij (10) 

Therefore Theorem 4 shows that if 



-< x ± q => p{M{x)) < 1 (11) 

then (9«r is the unique carrying simplex and T\F is a diffeomorphism. From (7), (8) and 
(10) we see that (11) holds in case one of the following conditions is satisfied: 

B ■ 

- 1 A iJ < 1 fOT alH > ( 12 ) 

A u . 

or 

^i^A, <l for allj (13) 
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These conditions thus imply a unique carrying simplex, by Theorem 4. 
To arrive at a biological interpretation of (12), we rewrite it as 

q^AjKl (14) 

j 

where Oj := -x*- is the axial equilibrium for species i, that is, its stable population in the 
absence of competitors. Equation (9) tells us that is the logarithmic rate by which the 
growth of population i inhibits the growth rate of population j. Thus (14) means that the 
average of these rates must be rather small compared to the single species equilibrium for 
population i. The plausibility of this xlis left to the reader, as is the biological meaning of 
(13). 

When n = 1, Equation (9) defines the map Tx = xe b ax of Example 2. The positive 
fixed point is q = a/b, and both (12) and (13) boil down to b < 1, which was shown 
to imply a unique carrying simplex. That example also showed that there is no carrying 
simplex when b > 2. As Equation (9) reduces to Example 2 on each coordinate axis, we see 
that Equation (9) lacks a carrying simplex provided 



- ^2 Aij > 2 for some i, 



A n j 



or 

^ A^ > 2 for some j 



Bi 



A. 



Example 6 

Consider a competing population model due to Leslie & Gower [19]: 

T:K^K, T iX = - , d, Aj > (15) 

1 + Z^j Ji ij x j 

Note that T need not be locally injective. When n = 1 all trajectories converge to if 
C < 1, and all nonconstant trajectories converge to if C > 1. The case n = 2 is 
thoroughly analyzed by Cushing et al. [6]. 
Here 

hence (C5) holds. We assume Ci > 1, guaranteeing (C4) with qi = ^ = ^-. In (5) we have 

1 + Ei A U x l 

so the row sums of M{x) are < 1 for all x provided % Ylj A ij < 1. Therefore when 



Mij(x) = — — — - < XiAij. 



1<Q<1 + ^ 



Theorems 3 and 4 yield the following conclusions: There is a global attractor T C [0, q], the 
unique carrying simplex is <9i<r, and T\T is a diffeomorphism. 
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Example 7 

Consider a recurrent, fully connected neural network of n cells (or "cell assemblies", Hebb 
[11]). At discrete times t = 0,1,..., cell % has activation level Xj(i) > and the state of 
the system is x(t) := (x\(t), . . . ,x n (t)). Cell i receives an input signal Si(x(t)) which is a 
weighted sum of all the activations plus a bias term. Its activation is multiplied by a positive 
transfer function n evaluated on Sj, resulting in the new activation Xi(t + 1) = Xj(i)Tj(sj). 

We assume each cell's activation tends to decrease the activations of all cells, but each 
cell receives a bias that tends to increase its activation. We model this with negative weights 
— Aij < 0, positive biases Bi > 0, and positive increasing transfer functions. For simplicity 
we assume all the transfer functions are e CT where a: [0, oo) — > [0, oo) is C 1 . States evolve 
according to the law 

T:K^K, Ti(x) = Xiexpa(si(x)), Sj(x) := Bj - y^AjjXj 

3 

We also assume 

cr(0)=0, a'{s)>0, supa'(s) = 7 < oo, (s € R) (16) 
It is easy to verify that (CI), (C2), (C4) and (C5) hold, with 

Qi ■= ~7~ j Gi(x) := expa(Bj - ^ AjXj), = a'(si(x))A i:j (17) 

where Mij(x) is defined as in (5). 

It turns out that for given weights and biases, the system has a unique carrying simplex 
provided the gain parameter 7 in (16) is not too large. It suffices to assume 

-1 

(18) 

For then (16), (17), (18) imply (8) and hence (C3), so Theorems 3 and 4 imply a unique 
carrying simplex for T. 

There is a vast literature on neural networks, going back to the seminal book of Hebb [11]. 
Network models of competition were analyzed in the pioneering works of Grossberg [7] and 
Cohen & Grossberg [5]. Generic convergence in certain types of competitive and cooperative 
networks is proved in Hirsch [15]. Levine's book [20] has mathematical treatments of several 
aspects of neural network dynamics. 

Competitive differential equations 

Consider a periodic differential equation in K: 

u i = UiGi(t,ui,...,u n ) = u i Gi(t + l,u 1 ,...,u n ), t,Ui>0, (i = l,...,n)) (19) 

where the maps Gi : K — > R are C 1 . The solution with initial value u(0) = x is denoted 
by t 1 ► T t x. Solutions are assumed to be defined for all t > 0. Each map T t maps K 



8 



diffeomorphically onto a relatively open set in K that contains the origin. The Poincare 
map is T :=T\. 

We postulate the following conditions for Equation (19): 

(Al) total competition: ^ i < 0, = 1, . . . , n) 

(A2) strong self-competition: J2kei(x) T^^^x) < 

(A3) decrease of large population: Gi(t,x) < for Xj sufficiently large. 
This implies existence of a global attractor for the Poincar'e map T. 

(A4) increase of small populations: Gi(t,0) > 0. 

Under these assumptions there are two obvious candidates for a carrying simplex for T, 
namely d^B and <9|<r, the respective boundaries in K of -B(O) and T. Existence of a unique 
carrying simplex implies d^B = d^T. 

Theorem 8 Assume system (19) has properties (Al) — (A4). Then there is a unique car- 
rying simplex, and it is the global attractor for the dynamics in K \ {0}. 

The proof, which will be given elsewhere, uses a subtle dynamical consequence of competi- 
tion discovered by Wang & Jiang [34, page 630]: If u(t),v(t) are solutions to Equation (19) 
such that for all i 

Ui(t) < Vi(t), (s <t < Si), 

then 

Example 9 

A competitive, periodic Volterra-Lotka system in K of the form 

Ui = Ui(Bi(t) - ^2Aij(t)ui), Bi,Aij > 
3 

satisfies (Al) — (A4) and thus the conclusion of Theorem 8. 
Example 10 

Several mathematicians have investigated carrying simplex dynamics for competitive, au- 
tonomous Volterra-Lotka systems in K having the form 

iij = Uj (Bj - AjjUj) := UjHj(ui, u n ), Bi,Aij>0 (20) 

3 

The best results are for n = 3: the interesting dynamics is on a 2-dimensional cell, therefore 
the Poincare-Bendixson theorem [10] precludes any kind of chaos and makes the dynamics 
easy to analyze. The dynamics for generic systems were classified by M.L. Zeeman [39], with 
computer graphics exhibited in Zeeman [40] . She proved that in many cases simple algebraic 
criteria on the coefficients determine the existence of limit cycles and Hopf bifurcations. 
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Van den Driessche and Zeeman [33] applied Zeeman's classification to model two com- 
peting species with species 1, but not species 2, susceptible to disease. They showed that if 
species 1 can drive species 2 to extinction in the absence of disease, then the introduction of 
disease can weaken species 1 sufficiently to permit stable or oscillatory coexistence of both 
species. 

Zeeman Sz Zeeman [42] showed that generically, but not in all cases, the carrying simplex 
is uniquely determined by the dynamics in the 2-dimensional facets of K. Systems with 
two and three limit cycles have been found by Hofbauer Sz So [18], Lu Sz Luo [21]), and 
Gyllenberg et al. [8]. No examples of Equation (20) with four limit cycles are known. 

More information on the dynamics of Equation (20) can be found in [32, 36, 37, 41, 43]. 

Background 

In an important paper on competitive maps, Smith [29] investigated C 2 diffeomorphisms 
T of K. Under assumptions similar to (CO) — (C5) he proved T is strictly retrotone and 
established the existence of the global attractor T and the repulsion basin B(0). He showed 
that <9k-B(0) and d^F are compact unordered invariant sets homeomorphic to the unit 
simplex, and each of them contains all periodic orbits except the origin. His conjecture that 
d^B = <9«r remains unproved from his hypotheses. He also showed that for certain types 
of competitive planar maps every bounded trajectory converges, extending earlier results of 
Hale Sz Somolinos [9], de Mottoni Sz Schiaffino [27]. 

Using Smith's results and those of Hess Sz Polacik [12], Wang Sz Jiang [34] obtained 
unique carrying simplices for competitive C 2 maps. 

For further results on the smoothness, geometry and dynamics of carrying simplices, see 
[1, 23, 24, 25, 26]. 

Mea culpa Uniqueness of the carrying simplex for Equation (20) was claimed in Hirsch 
[13], but M.L. Zeeman [38] discovered an error in the proof of Proposition 2.3(d). 
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